Numerical Renormalizat ion- Group Study of Particle-Hole Symmetry Breaking 

in Two-Channel Kondo Problem : 
Effect of Repulsion between Conduction Electrons and Potential Scattering 
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Particle-hole symmetry breaking perturbation in two-channel pseudospin Kondo problem is studied 
by the numerical renormalization-group method. It is shown that the repulsion between conduction 
electrons at the impurity site and the single particle potential are the relevant perturbations against 
the conventional non- Fermi liquid fixed point. Although the repulsion (potential) with realistic 
strength prevents the overscreening of pseudospin, it induces in turn a real spin, which is also over- 
screened again. Thus the real spin susceptibility becomes anomalous contrary to the conventional 
two-channel Kondo problem. 
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Multichannel Kondo problem has attracted much at- 
tentions in these years because of its anomalous non- 
Fermi liquid behavior. While the problem was origi- 
nally discussed long ago as a generalized Kondo effect 
of magnetic impurity with orbital degeneracy Q, the 
two-channel Kondo problem has recently been revived 
in a proposal of quadrupolar Kondo effect as an origin of 
U-based heavy fermions The two-level system inter- 
acting with conduction electrons was also recognized as a 
candidate for realization of two-channel Kondo model || . 
The latter system has attracted much interest not only 
because it offers a model explaining anomalous transport 
properties of glassy metals (||| but also it is expected 
to give a canonical model of strong coupling electron- 
phonon systems 

Although the two-channel Kondo problem has been 
fully solved by a variety of methods |^-|l5|, there seems 
still to remain for us to clarify a reality of the model it- 
self. As for relation with experiments, the pseudospin 
two-channel Kondo model gives anomalous behavior the 
susceptibility of the pseudospin (i.e., that of charge po- 
larization), while it is the real spin susceptibility that 
exhibits non-Fermi liquid behaviors in some of the above 
systems. The latter point may be related to an appear- 
ance of localized real spin. For instance, a repulsion be- 
tween conduction electrons at the impurity site, which 
was neglected in the above pseudospin models |ll| , is ex- 
pected to prevent the overscreening p7[ and to induce a 
real spin. 

The purpose of this paper is to examine the effect of 
such a repulsion on the two-channel Kondo model by nu- 
merical renormalization-group (NRG) method (l§[LSj. It 
is shown that the fixed point Hamiltonian is described 
not only by the conventional exchange coupling J* but 
also an impurity potential V* , to which the repulsive in- 
teraction U is renormalized. Namely, the single particle 
potential V is also a relevant perturbation. The compe- 



tition between the exchange coupling and the repulsion 
or the impurity potential induces the degrees of freedom 
of channel (i.e., real spin) and leads to the pseudospin 
singlet ground state for realistic strength of U or V. It 
is the particle-hole symmetry breaking that causes such 
competition. The overscreening of so induced real spin 
makes again the real spin susceptibility anomalous con- 
trary to the conventional pseudospin two-channel Kondo 
problem without U and V. 

We begin with the model Hamiltonian for the Wilson 
NRG calculation as follows: 
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where the indices m and a denote a label indicating chan- 
nel 1, 2 and pseudospin f, J., respectively, er is the Pauli 
matrix vector, and t is twice the operator of the impurity 
pseudospin of 1/2. The pseudospin stands for the charge 
polarization at the impurity site and two channels for two 
components of "real" spin of conduction electron. Here 
we have defined 
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where 2D denotes the bandwidth of conduction electrons, 
J the exchange interaction between conduction electrons 
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and impurity pseudospin, V the potential at impurity 
site, and U the repulsion between conduction electrons 
at impurity site Hereafter we set D = 1, i.e., the 

energy levels are scaled by D and ignore A-dependence 
in £„, i.e., £ n = 1, because £ n — > 1 for large n. 

The conserved quantities of the Hamiltonian Hn, (Q), 
are the total number of conduction electrons Q, the real 
spin j and the total pseudospin S, defined as follows: 
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Since both the repulsion U and the potential V breaks 
the particle-hole symmetry unless 3t//2 + V — 0, the de- 
generate eigenstates denoted by ±Q are split in general. 
In our calculations, we have used A = 3 and retained low 
lying energy states up to 300 states at each step as bases 
for constructing new quadruple states. 

First we have investigated the case U = V = and 
verified that the same energy levels as the work of Pang 
and Cox jl0| are reproduced. Next we have investigated 
the case U ^ 0. The flow diagram of level of low lying 
states for J = 2.0 and U = 1.6 is shown in Fig. [TJ. The 
solid (dotted) lines are for even (odd) iterations. Each 
level is labelled by (Q, j, S). The ground state of the fixed 
point is pseudospin doublet (S = 1/2), which is expected 
for the case where the exchange coupling is stronger than 
the repulsion. In Fig. |^ the flow diagram for J = 1.0 and 
U = 2.0 is shown. The ground state is now pseudospin 
singlet (S — 0), because the repulsion U, larger than 
exchange coupling J and the hopping D = 1, prohibit the 
overscreening. It is noted that the ground state is still 
degenerate due to degrees of freedom of channel, i.e., j = 
1/2. It is remarked that the position of energy levels at 
the fixed point in Fig. [I] and Fig. || exactly coincide with 
each other while the nature of ground state is different. 
As discussed later, the quantum numbers specifying each 
levels in Fig. ^ are in a certain relation with those in Fig. 

The nature of ground states for various coupling con- 
stants J, U and V = are shown in Fig. |^. The closed 
circles stand for the ground state with S = 0, while the 
open circles S = 1/2. The line dividing two types of 
ground state is drawn by estimating the coupling con- 
stants which give the same energies of these two types 
of ground state. It is noted that the boundary line flat- 
ten as J — > and has a constant slope for J>1. We 
can understand this result as follows. The energy gains 



for overscreening formation at the impurity (n = 0) site 
are both due to the exchange energy J and the kinetic 
energy associated with the transfer D between 0- and 1- 
site, while the energy loss arises through the repulsion 
U between overscreened conduction electrons. On the 
other hand, for singlet formation, the energy gain is due 
to avoiding the repulsion U while the energy loss is due 
to the exchange coupling J and the kinetic energy D. 
Consequently, the boundary line is roughly determined 
by the condition U ~ max(J, D). Namely, for U > J/8 
and D/4, the ground state becomes a pseudospin singlet. 
It is noted that the ground state is expected to belong to 
that of S = for a realistic value of U and J. 

Now we discuss about properties of the fixed point. 
The fixed point Hamiltonian H* is described as 

oo 
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where H* nt has the same symmetry as -ffi OC ai, in (J]]), the 
effective couplings J* , a*, /?*, and V* may depend on 
the initial couplings in general, and e is a constant energy 
shift. If we set J* = J, a* = f3* = -U, V* = 3C//2 + V, 
and e = 3U + 2V, H* nt becomes equivalent to H\ oca \ in 

Since the energy of low lying excited states at the fixed 
point are mainly determined by H* nt , we can determine 
the parameters J* ~ e in (0) so as to reproduce the low 
lying energy levels at N = 39. The results for the initial 
parameters, (J,U,V), (a) (0.5,0.0,0.0), (b) (2.0,1.6,0.0), 
c) (1.0,2.0,0.0) and (d) (0.2,0.4,-0.6) are shown in Table 
It is noted that the effective exchange coupling, J* is 
independent of the initial coupling U (including U = 0), 
and a* and (3* are always zero. In the case (d), there is 
the particle-hole symmetry so that the fixed point is the 
same as the case (a) where U = V = 0. The character 
of the fixed point is determined mainly by the effective 
impurity potential V* which depends on the initial cou- 
plings J, U and V, i.e. V* — f(J,U,V). Consequently, 
the effective interaction at the fixed point can be written 
as 
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(J* = 0.20). (8) 



J, [/-dependence of V* with V = is shown in Fig. [|. It 
is noted that V* increases (decreases) as U (J) increases. 
From this effective interaction, the "flow lines" for scaling 
in parameter space are obtained from /( J, U, V = 0) = 
constant. Especially, for V* = J* = 0.20, the "flow line" 
becomes equivalent to the boundary line shown in Fig. 
|3|, because the first excitation energy is zero for these 
couplings. 

In order to investigate to and T dependence of suscep- 
tibility, let (Q,j, S) be (QdJd, S d ) for V = V£< 0.20, 
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where the ground state is pseudospin doublet (S = 1/2), 
and (Qs,js,Ss) for V* = Vg > 0.20 where the ground 
state is pseudospin singlet (S = 0). For each low lying 
excited states, we can find the relations, Qg = —Qd — 1, 
is = S D , and Sg = 3d- 
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the low lying excited energies at the fixed point for each 
parameters are the same as easily seen by means of the 
effective interaction (^). A prime example is the relation 
between Fig. |l| and Fig. || as mentioned above. Ac- 
cording to this example, the coincidence of energy levels 
occurs after 20 iterations. From this coincidence, it is ex- 
pected that u and T(uj/D,T/D < A~ 20 / 2 ) dependence 
of the susceptibility of real spin (channel) for V£ coin- 
cide with those of the pseudospin for Vg — 2 J* — V£ and 
vice versa. This is a new aspect of two-channel Kondo 
problem which has not been recognized so long as the 
conventional model without the repulsion and the poten- 
tial scattering (U = V = 0) had been investigated. 

This remarkable aspect can be seen more vividly by 
investigating the spectral weight of the dynamical sus- 
ceptibilities for real spin of conduction electrons at the 
impurity site, ancl for the impurity pseudospin, 

X"(w). They are calculated by the method of ref. as 
shown in Figs. || and |[ It is noted that, in the presence 
of the repulsion U, Xj( w ) shows the non- Fermi liquid be- 
havior with lim^^o ( w ) being finite, while without the 
repulsion it shows Fermi liquid one with x'J (0) = 0. How- 
ever, if we set 3U/2+V = 0, the similar calculations show 
that x'j (0) = 0. Namely, it is the particle-hole symmetry 
breaking that gives x"(w)'s the non- Fermi liquid behav- 
iors. The potential V* shifts the number of conduction 
electrons at the impurity site from one in each channel, 
though the exchange works to hold the overscreening for- 
mation. This competition induces the degrees of freedom 
of channel (i.e., real spin). It is also overscreened again by 
conduction electrons with two channels, i.e., pseudospin 
degrees of freedom. Thus the real spin susceptibility be- 
comes anomalous due to the potential V* which breaks 
the particle-hole symmetry. The case V* < is under- 
stood as V* > by particle-hole transformation. Since 
the degrees of freedom of pseudospin, however, are not 
perfectly vanished, pseudospin susceptibility is anoma- 
lous for any strength of the repulsion including U = 0. 

In summary, the low lying excited states at the fixed 
point of pseudospin two-channel Kondo model with the 
particle-hole symmetry breaking perturbations are de- 
scribed not only by the exchange J* but also the impurity 
potential V* . For |V~* > 0.20, i.e., U > J/8 and D/4, 
realistic values, a pseudospin singlet ground state is real- 
ized in contrast with the pseudospin doublet ground state 
which is realized in the conventional two-channel Kondo 
problem. The spectral weight of the dynamical suscepti- 
bility of the real spin shows the non- Fermi liquid behav- 
ior because of the overscreening of the real spin. Thus, 
it is expected that the magnetic non-Fermi liquid behav- 
iors observed in some compounds are understood by the 



particle-hole symmetry breaking perturbation, which in- 
duces the degrees of freedom of the real spin. 
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an isotropic one by renormalization-evolution as will be 
discussed elsewhere; so that we confine ourselves to dis- 
cussing the case of isotropic repulsive interaction in (2). 

FIG. 1. The flow diagram for J = 2.0, U = 1.6 and V = 0. 
Each level is labelled by (Q,j,S), where Q, j, and S are 
the total number of electrons, the real spin, and the total 
pseudospin, respectively. The solid (dotted) lines are for even 
(odd) iterations. The ground state is a pseudospin doublet. 

FIG. 2. The flow diagram for J = 1.0, U = 2.0 and V = 0. 
The ground state is a pseudospin singlet. The degeneracy of 
the ground state is due to degrees of freedom of channel, i.e., 
3 = 1/2. 

FIG. 3. The nature of ground states for various sets of 
coupling constants of J, U and V = in the unit D. The 
open circles are for the ground state with pseudospin doublet 
S = 1/2, while the closed circles the ground state with pseu- 
dospin singlet S = 0. The line dividing two phases of, S = 
and S = 1/2, is drawn by inspection. 

FIG. 4. J, (7-dependence of the effective impurity poten- 
tial, V* , with V = at the fixed point. 

FIG. 5. ui dependence of imaginary part of local real spin 
susceptibility x'j 

FIG. 6. io dependence of imaginary part of impurity pseu- 
dospin susceptibility \t (u>). 





(J,U,V) 


J* 


Q* 




V* 


e 


(a) 


(0.5,0.0,0.0) 


0.20 











0.80 


(b) 


(2.0,1.6,0.0) 


0.20 








0.12 


0.80 


(c) 


(1.0,2.0,0.0) 


0.20 








0.28 


0.88 


(d) 


(0.2,0.4,-0.6) 


0.20 











0.80 



TABLE I. Effective couplings, J* ~ e, which make a re- 
production of the energy levels at iV = 39 for the initial 
parameters, (J,U,V), (a) (0.5,0.0,0.0), (b) (2.0,1.6,0.0), (c) 
(1.0,2.0,0.0), and (d) (0.2,0.4,-0.6). 
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